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Abstract 

We have calculated the single jet inclusive cross section as measured at Fermilab in 
next-to-leading order QCD using recent parton distributions of the CTEQ collaboration. 
We studied the scheme dependence of the jet cross section by employing the MS and DIS 
factorization schemes consistently. For Et > 200 GeV, we find that the cross section in 
the DIS scheme is larger than in the MS scheme yielding a satisfactory description of the 
CDF data over the whole Et range in the DIS scheme. 



Recently, the CDF collaboration presented a precise measurement of the inclusive differential 
cross section for jet production in pp collisions at 1.8 TeV [3J. The measurement was compared 
to next-to-leading order (NLO) perturbative QCD predictions @ for jet transverse energies Et 
in the range of 15 to 440 GeV in the central pseudorapidity region 0.1 < \r]\ < 0.7, using a 
selection of commonly used parton distributions. The experimental results for Et > 200 GeV 
show evidence of a possible deviation as compared to the NLO prediction based on the current 
sets of parton distributions, which are obtained from global analyses || f|, || || of deep inelas- 
tic lepton-nucleon scattering and related data. Before explanations for this deviation based on 
new physics [[?[ are seriously considered, it is crucial to study possible explanations within the 
Standard Model. 



Clearly, the perturbative predictions of the jet cross section depend on the parton distri- 
butions. Then the question arises whether the parton distributions can be adjusted to accom- 
modate the jet measurements while still having a good description of the data sets as used in 
previous global analyses. This vital question has been addressed by two groups. Huston et al. 
H have carried out a new global QCD analysis incorporating the CDF inclusive jet data with 
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E T > 75 GeV (the very precise data in the low E T range are excluded due to potential theo- 
retical and experimental problems) and reached the conclusion, that there is enough flexibility 
in NLO global analyses to enhance the high inclusive jet cross section by 25 — 35% above 
previous calculations. A similar analysis was done by Glover et al. ||. They also incorporated 
the CDF single jet inclusive measurement in a global NLO parton analysis of the available deep 
inelastic and related data. In particular, they included the jet data in the region Et > 50 GeV 
in their fit. They find that it is impossible to accommodate both the jet data over the complete 
Et range and the deep inelastic structure function data. However, the CDF data for Et < 200 
GeV and the deep inelastic data were reasonably well fitted simultaneously. A fit over the full 
Et range yields quarks which are completely incompatible with the large-x structure function 
data, a value for a s , which is larger than the current world average, and a renormalization of 
the CDF jet data, which is barely compatible with the allowed range So from these two 
studies it still seems to be unclear, whether the difference between the CDF inclusive jet cross 
section data and the NLO predictions can be attributed to unsufficient knowledge of the parton 
distributions. 



It is well known that there are other uncertainties of the perturbative calculations which 
should be considered. First, we have the dependence on the renormalization and factorization 
scales. Second, changes in the strong coupling a s resulting in changes of Aqqq must be con- 
sidered. These uncertainties have been discussed in || and mainly affect the renormalization 
and not the shape of the inclusive jet cross section. Third, there exists the dependence of the 
perturbative predictions on the factorization scheme. To our knowledge, this dependence of 
the single jet cross section has not been investigated yet. 



Beyond leading order, the parton distribution functions f a (x, fi 2 ) depend on the factorization 
scheme. In applications to physical processes, the scheme chosen for the parton distributions 
must match that for the hard scattering cross section in the parton model formula The 



same parton distributions in two different schemes differ because of two effects. First, the 
scheme influences the evolution kernel. Second, the functional dependence at the starting scale 
(i = Qo differs since the parton distributions are fitted to the same data with the evolution 
kernel for the appropriate scheme. Thus, they are functionally equivalent in the sense that they 
yield the same physical cross sections for the data included in the analysis. The MS parton 
distributions are guaranteed to satisfy the momentum sum rule. The DIS scheme [|TTJ was 
defined specifically to make the relation between the parton distributions and F% as simple as 
possible by absorbing all the NLO terms into the definition of f% IS ( x i A* 2 )- This does not yet 
define the DIS prescription for the gluon distribution. Here, the convention has beed adopted 
to require the momentum sum rule to be preserved in the DIS scheme as well. However, this 
requirement fixes only the second moment of the gluon distribution, so that one completes the 
definition by requiring the condition on the second moment to be valid for all moments. This 



convention was first introduced by Diemoz et al. |12[ and is now generally applied for parton 
distribution function analyses in the DIS scheme. 



Functionally equivalent parton distributions in the MS and the DIS schemes have been 
constructed by the CTEQ collaboration |4], [5]. The most recent ones are the sets CTEQ3M 
(MS) and CTEQ3D (DIS), respectively ||. These two sets are obtained from independent fits 
to the same data sets under the same assumptions except for calculating the evolution kernel. 
The total \ 2 f°r the two sets is supposed to be very similar. Unfortunately, this is not stated 
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explicitly in ||. However, this fact was reported for the earlier sets CTEQ1M and CTEQ1D 
[|J. Thus, the two sets describe totally equivalent physics as defined by the data sets selected 
in the analysis. Now we can ask ourselves whether they also yield the same single jet inclusive 
cross sections or not. 



There also exist DIS parton distributions from the other collaborations || . In particular, 
MRS(DO') and MRS(D-') sets in the DIS scheme are available H|. They are obtained from the 
corresponding sets in MS, which were fitted to the deep inelastic and other data, by applying 
the 0(a s ) perturbative transformation formula between the two schemes. It is known that this 
may be unreliable in situations where the NLO terms involving, for example, a large gluon 
contribution, are of comparable size as the LO term involving small sea quarks [||. Further- 
more, the DIS version of MRS is not equivalent to the Diemoz et al. description, since only 
the quark distribution functions are transformed. The GRV collaboration |J also constructed 
a DIS set in their recent 1994 analysis (GRV(94)) by fitting MS and DIS parton distributions 
independently to the deep inelastic and other data. Unfortunately, this analysis is restricted 
to three flavors with the charm contribution to the deep inelastic structure function generated 
perturbatively. This makes this set less useful for calculations of other hard scattering cross 
sections. Therefore, we shall restrict ourselves to the CTEQ3 sets in the following. 



Our calculation uses the NLO parton level program JETSAM which originally was 
designed for the calculation of inclusive single jet cross sections in resolved photoproduction. 
This program was developed on the basis of the theory of Aversa et al. ]15[] which uses the 
phase space slicing method to cancel collinear divergences. The cuts and the algorithm for 
defining jets out of up to three partons are modelled as closely as possible to the experimental 
set-up. The jets are defined according to the Snowmass algorithm jTJj] with a jet cone size of 
R = 0.7 and lie in the pseudorapidity range between 0.1 < \r]\ < 0.7 as in the CDF analysis. 
The factorization and renormalization scales are chosen to be fip = fiR = -Et/2. The Monte 
Carlo JETSAM contains the necessary modifications of the NLO hard scattering cross sections 
to run it in the MS and the DIS scheme as defined above. 



First, we compare the CDF jet data with the NLO prediction obtained from the CTEQ3M 
set of partons. The fractional difference between the data and the theoretical prediction is 
plotted in Fig. 1. The data are taken from Table 2 in 0. Over the observed range of Et, the 
experimental cross section decreases by more than a factor 10 8 and the quadratic sum of the 
correlated systematic uncertainties grows from ±18% at Et = 50 GeV to ±36% at Et = 400 
GeV. In Fig. 1, the agreement of the data with the theory for the CTEQ3M set is measured 
by the distance of the data points from the horizontal line at zero. The normalization shown 
is absolute. These results show excellent agreement in shape and normalization for Et < 200 
GeV, where the cross section falls by six orders of magnitude. Above Et = 200 GeV, the 
CDF cross section is significantly higher than the NLO cross section. This agrees with the 
statements made in ref. [fl], Bj, where other parton sets in the MS scheme obtained from 
recent global analyses had been selected for the NLO calculation of the jet cross section. Thus, 
the agreement below 200 GeV and disagreement above 200 GeV seems to be a unique feature 
of all MS parton sets from recent global analyses which include both fixed-target and HERA 
deep inelastic data. 
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Figure 1: Fractional differences of the CDF inclusive jet data, the CTEQ3D prediction, and the 
CTEQ3D prediction without compensation terms from the CTEQ3M NLO QCD predictions. 



In Fig. 1, we have also plotted the prediction for the CTEQ3D parton distributions. It is 
given by the dashed line which presents the fractional difference of the CTEQ3D cross section 
from the CTEQ3M cross section, i.e. (CTEQ3D-CTEQ3M) /CTEQ3M. As we can see, this 
line deviates very little from the horizontal line (the CTEQ3M result) for E? < 150 GeV, 
showing that the single inclusive jet cross section does not depend on the chosen factorization 
scheme in this Et range. Above Et = 150 GeV, the DIS prediction starts to deviate from the 
horizontal line and increases up to 0.4 at Et = 400 GeV. Thus, in the large Et range the jet 
cross section depends on the factorization scheme. For Et > 150 GeV, the DIS cross section is 
larger than the MS cross section so that it agrees with all the large Et data points except the 
last one which, however, has a large statistical error. In addition, we must take into account 
the systematic error of the experimental data Pj] . 

In Fig. 1, we also show a second curve (dotted). This presents the fractional difference to 
the CTEQ3M prediction, obtained also with the DIS parton distribution set CTEQ3D, but 
without the appropriate modification of the hard scattering cross section to the DIS scheme, 
i.e. this cross section is left as calculated in the MS scheme. Of course, this is an inconsistent 
procedure. But it shows that the deviation of the DIS cross section from the MS cross section 
at Et = 400 GeV, which is 40%, comes to a large extent from the modified parton distributions 
(25%) and to a lesser extent from the modified subtraction terms in the hard scattering cross 
section (15%). At the larger values of Et, these two effects do not compensate each other as 
one might expect. At the smaller Et < 150 GeV, this compensation occurs, the dashed curve 
is nearer to the horizontal line than the dotted one, although the difference between the two 
curves is very small. 
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We also performed the same calculations for the older sets CTEQ2M and CTEQ2D ||, 
respectively. The results are similar. The agreement between the data and the predictions 
with either the MS or DIS parton distributions in the Et < 150 GeV is less satisfactory. The 
fractional difference of the CTEQ2D cross section from the CTEQ2M cross section shows the 
same increase with increasing Et up to the value 0.5 at Et = 400 GeV. 

Fig. 2 compares the quark and gluon distributions at Q 2 = 10 and 10 4 GeV 2 of the two sets 
CTEQ3M and CTEQ3D. We plot the singlet quark distribution J2 a (Qa + Q a ) and the gluon dis- 
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Figure 2: The CTEQ3M and CTEQ3D parton distributions J^aiQa + Q a ) an d 9 a t Q 2 = 10 and 
10 4 GeV 2 . 

tribution g multiplied by x 2 , so that the area under the curves is the total momentum fraction 
carried by the quarks and gluons. We see that for both Q 2 values, the DIS quarks are larger 





5 



(smaller) than the MS quarks for x > 0.3(x < 0.3), whereas the DIS gluons are larger than the 
MS gluons over the whole x range of 0.01 < x < 1. At smaller Et, where the jet cross section 
receives contributions mainly from qg— and gg— initiated subprocesses |J, the compensation 
between the change of the parton distributions from MS to DIS and the modification of the 
hard scattering subprocesses seems to work reasonably well (see Fig. 1). For Et > 200 GeV, the 
jet production is dominated by qq initiated subprocesses ||. Since centrally produced jets of 
transverse energy Et sample partons at x ~ xt = ^t 21 , we can understand that the increased 
quark distribution for x > 0.3 in the DIS scheme leads to an increased jet cross section for 
Et > 250 GeV in case that in this region the compensation is distorted due to the singling 
out of a particular subprocess. The increase of the gluon distribution in this region is small, so 
that the qq process is mainly responsible for the increase of the jet cross section. For smaller 
Et, the effect of the larger gluon is compensated by a smaller quark distribution for x < 0.3. 

In summary, we have studied the dependence of the inclusive jet cross section on the factor- 
ization scheme. We have considered two schemes, MS and DIS, for which parton distributions 
have been constructed in a global analysis by the CTEQ collaboration. We find that the inclu- 
sive jet cross section for Et > 200 GeV in the DIS scheme is up to 40% larger than in the MS 
scheme and accounts well for the CDF inclusive jet cross section data in this Et region as well 
as in the low Et range. 

Of course, the significant scheme dependence of the high Et jet cross section is an artifact of 
finite order perturbation theory. In infinite order, this dependence is absent. From this point of 
view, our results indicate that the NLO theory is inadequate and one must go to the NNLO in 
order at least to diminish the scheme dependence and to obtain more reliable results. However, 
since parton distribution functions are involved, which are determined from other processes, 
the problem of the scheme dependence might have another solution. It is conceivable, that the 
parton distributions in the DIS scheme are the right ones with enough flexibility built in to 
describe all deep inelastic data together with the CDF jet data. Then, the task would be to 
find the corresponding MS set which describes the same data by perhaps changing the input 
parametrization of the distribution functions, in particular that of the gluon, which is much 
less constrained by deep inelastic data than the quark distribution. Such an investigation is 
beyond the intention of this work. We therefore conclude that it might still be possible that 
the difference between the CDF inclusive data and the NLO QCD predictions can be due to a 
deficiency in the way parton distributions in the MS scheme are constructed in the usual global 
analyses work as performed by the CTEQ, MRS, and GRV collaborations. 
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